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I. INTRODUCTION 



In the leading logarithmic approximation, the high-energy scattering in perturbative 
QCD is determined by the BFKL pomeron [|TJ. It is well known that the power behavior 
of BFKL cross section violates the Froissart bound. The BFKL pomeron describes only 
the pre-asymptotic behavior at not very large energies and in order to find the true high- 
energy asymptotics in perturbative QCD we need to unitarize the BFKL pomeron. This is 
a difficult problem which has been in a need of a solution for more than 20 years. However, 
until recently, it was a common belief that at least at present energies (e.g. for small-x 
deep inelastic scattering in HERA) the corrections to BFKL pomeron are small so they 
can be neglected. Contrary to that expectations, recent calculation of the next-to-leading 
correction to the BFKL kernel || shows that this correction is very big. It is very likely 
that further corrections are also large which means that we must deal with the problem of 
the unitarization of the BFKL pomeron at present energies. 

One of the most popular ideas on solving this problem is to reduce the QCD at high 
energies to some sort of two-dimensional effective theory which will be simpler than the 
original QCD, maybe even to the extent of exact solvability. Some attempts in this direction 
were made starting from the work [|J but the matter is an open issue for the time being. 
In this paper I will describe the new approach to the effective action which is based on the 
factorization in rapidity space for high-energy scattering. 

The form of factorization is dictated by process kinematics (for a review, see ||). A 
classical example is the factorization of the structure functions of deep inelastic scattering 
into coefficient functions and parton densities. In the case of deep inelastic scattering, there 
are two different scales of transverse momentum and it is therefore natural to factorize the 
amplitude in the product of contributions of hard and soft parts coming from the regions of 
small and large transverse momenta, respectively. On the contrary, in the case of high-energy 
(Regge-type) processes, all the transverse momenta are of the same order of magnitude, but 
colliding particles strongly differ in rapidity so it is natural to factorize in the rapidity space. 
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Factorization in rapidity space means that the high-energy scattering amplitude can be 
represented as a convolution of contributions due to "fast" and "slow" fields. To be precise, 
we choose a certain rapidity 770 to be a "rapidity divide" and we call fields with 77 > r] fast 
and fields with rj < i]q slow where i]q lies in the region between spectator rapidity tja and 
target rapidity rjB- (The interpretation of this fields as fast and slow is literally true only in 
the rest frame of the target but we will use this terminology for any frame). 

To explain what we mean by the factorization in rapidity space let us recall the operator 
expansion for high-energy scattering || where the explicit integration over fast fields gives 
the coefficient functions for the Wilson-line operators representing the integrals over slow 
fields. For a 2=^2 particle scattering in Regge limit s ^> m 2 (where m is a common mass 
scale for all other momenta in the problem t ~ p\ ~ (p'a) 2 ~ Pb ~ (p'b) 2 ~ m2 ) we have: 



(As usual, s = (pa+Pb) 2 and t = (pa—p'a) 2 )- Here X{ (i = 1, 2) are the transverse coordinates 
(orthogonal to both pa and pb) and U{(x) = W (x)^-£^-U(x) where the Wilson-line operator 
U(x) is the gauge link ordered along the infinite straight line corresponding to the "rapidity 
divide" t]q. Both coefficient functions and matrix elements in Eq. (P depend on the 770 
but this dependence is canceled in the physical amplitude just as the scale \i (separating 
coefficient functions and matrix elements) disappears from the final results for structure 
functions in case of usual factorization. Typically, we have the factors ~ (g 2 In s/m 2 — i] ) 
coming from the "fast" integral and the factors ~ g 2 r] coming from the "slow" integral so 
they combine in a usual log factor g 2 In s/ m 2 . In the leading log approximation these factors 
sum up into the BFKL pomeron (for a review see ref. J7[). 

Unlike usual factorization, the expansion ([!]) does not have the additional meaning of 
perturbative v s nonperturbative separation - both the coefficient functions and the matrix 
elements have perturbative and non-perturbative parts. This happens due to the fact that 
the coupling constant in a scattering process is is determined by the scale of transverse 



A(Pa,Pb ^PaiP'b) 



(1) 



Efd 2 



x l ...d 2 x n C h - in (x 1 ,...x n )(p B \Tr{U il (x 1 )...U in (x n )}\p' B ) 
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momenta. When we perform the usual factorization in hard (k± > //) and soft (k± < 
fi) momenta, we calculate the coefficient functions perturbatively (because a s {k^_ > fi) is 
small) whereas the matrix elements are non-perturbative. Conversely, when we factorize the 
amplitude in rapidity, both fast and slow parts have contributions coming from the regions 
of large and small k±. In this sense, coefficient functions and matrix elements enter the 
expansion ([!]) on equal footing. We could have integrated first over slow fields (having the 
rapidities close to that of P b) and the expansion would have the form: 

=E/^z 1 ...d 2 z n I>-^^ (2) 

In this case, the coefficient functions D are the results of integration over slow fields ant the 
matrix elements of the U operators contain only the large rapidities rj > tjq. The symmetry 
between Eqs. (1) and (2) calls for a factorization formula which would have this symmetry 
between slow and fast fields in explicit form. 

I will demonstrate that one can combine the operator expansions ([I]) and (|2p in the 
following way: 

A(s,t) = ^2^J ^X^.^Xn (3) 
(p^|^ i i(x 1 )...^»(xOb^)(p B |^(x 1 )...C/f n "(x n )b , B ) 

where Uf = Tr(A a £7j) (A a are the Gell-Mann matrices). It is possible to rewrite this factor- 
ization formula in a more visual form if we agree that operators U act only on states B and 
B' and introduce the notation V% for the same operator as Ui only acting on the A and A' 
states: 

A(s,t) = (p A \(p B \exp(ifd 2 xV ai (x)U«{x)) \p' A )\p' B ) (4) 

In a sense, this formula amounts to writing the coefficient functions in Eq. ([!]) (or Eq. (|2|)) 
as matrix elements of Wilson-line operators. Eq. (Q) illustrated in Fig.l is our main tool for 
factorizing in rapidity space. 

In order to define an effective action for a given interval in rapidity T]q > r\ > r]' we use 
the master formula (0) two times as illustrated in Fig. 0. We obtain then 



FIG. 1. Structure of the factorization formula. Dashed, solid, and wavy lines denote photons, 
quarks, and gluons, respectively. Wilson-line operators are denoted by dotted lines and the vector 
n gives the direction of the "rapidity divide" between fast and slow fields. 

A(s,t) = (pA\(pB\e iS °«M> Y ^p' A )\p' B ) (5) 

where the Wilson-line operators Yi(x±) have the same form as Ui(x±) but aligned along the 
n' direction (and act only on B and B' states, cf. eq. (||)). In this formula, the region 
of rapidities greater than rj is represented operators Vi acting on the spectator A and A' 
states, the region of rapidities lower than rj' by the operators Y{ acting on target B and 
B' states, and the region 770 > rj > rj' Q is integrated out -all the information about it is 
contained in the effective action S e ff(Vi, Yj). As we shall see below, this effective action is in 
general non-local (unlike the local interaction term / dx±Vi(x±)Ui(x±) in the factorization 
formula (|])). Moreover, it contains the factors g 2 (r)o — rj' Q ) which are the usual high-energy 
logarithms g 2 ln^ where the energies so and s' correspond to rapidities 770 and 1]' . If we 
had a complete expression for S e g(rj ,rj' Q ) we could take i] = t/a (rapidity of the spectator 
particle) and r]' Q = t]b (rapidity of the target particle), then all the logarithmic dependence 
on the energy would be included in the effective action and the resulting matrix elements 
of the operators Vi between A states and operators Y between B states will contain no 
logarithms (and may me calculated in the first order in perturbation theory for a suitable 
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FIG. 2. The effective action for the interval of rapidities t]q > r] > The two vectors n and n' 
correspond to "rapidity divides" t]q and r]' Q bordering our chosen region of rapidities 

A and B particles such as virtual photons). Since multiple rescatterings are taken into 
account by e eff automatically the corresponding amplitude must be unitary. This program 
is probable not less difficult than the direct calculation of the many-pomeron exchanges in 
the perturbation theory but for the case of effective-action language we have some additional 
powerful methods such as semiclassical approach. 

The paper is organized as follows. In Sect. 2 we remind the Wilson-line operator language 
for small-x physics. The factorization formula is derived in Sect. 3 and in Sect. 4 we use 
it to define formally the high-energy effective action for a given interval in rapidity (Some of 
the results of this Sections were reported earlier in the letter ||). A semiclassical approach 
to calculation of this effective action is discussed in Sect. 5 and Sect. 6 contains conclusions 
and outlook. 

II. OPERATOR EXPANSION FOR HIGH-ENERGY SCATTERING 

Let us now briefly remind how to obtain the operator expansion For simplicity, 
consider the classical example of high-energy scattering of virtual photons with virtualities 

6 



~ — m 2 . 



A(s,t) = -i(0\T{j(p A )j(p' A )j(p B )j(p' B )}\0). (6) 

where j(p) is the Fourier transform of electromagnetic current multiplied by some suit- 
able polarization e M (p). At high energies it is convenient to use the Sudakov decomposition: 

p" = c^ + A^+K (7) 

where pi and p% are the light-like vectors close to pa and p B , respectively (p A = pi — 
P2Pa/ s > Pb = P2~ PiP 2 b/ s )- We want to integrate over the fields with a > o where o is 
defined in such a way that the corresponding rapidity is r] . (In explicit form r] = In | 
where a = — ). The result of the integration will be given by Green functions of the fast 
particles in slow "external" fields |J (see also ref. [Kj). Since the fast particle moves along 
a straight-line classical trajectory the propagator is proportional to the straight-line ordered 



gauge factor U [11]. For example, when x + > 0, y + < it has the form 



G(l , tf)=j y (WW5 ^. I7( ,J_^_ (8) 

We use the notations z, = z^pi and z* = z^p 1 ^ which are essentially identical to the light- 
front coordinates Z+ — Zst/y/s, Z— — Z m j\fs- The Wilson-line operator U is defined as 

U (x±) = [oopi + -oopi + (9) 

where [x, y] is the straight-line ordered gauge link suspended between the points x and y: 

[x,y]=Pexp (igfv du(x - y^A^ux + (1 - u)y)) (10) 

The origin of Eq. (|j) is more clear in the rest frame of the "A" photon (see Fig.2). Then 
the quark is slow and the external fields are approaching this quark at high speed. Due to 
the Lorentz contraction, these fields are squeezed in a shock wave located at = (in a 
suitable gauge like the Feynman one). Therefore, the propagator (|S[) of the quark in this 
shock-wave background is a product of three factors which reflect (i) free propagation from 



FIG. 3. Quark propagator in a shock-wave background 

x to the shock wave (ii) instantaneous interaction with the shock wave which is described 
by the operator U(z±) [] , and (hi) free propagation from the point of interaction z to the 
final destination y. 

The propagation of the quark-antiquark pair in the shock-wave background is described 
by the product of two propagators of Eq. (|8|) type which contain two Wilson-line factors 
U(z)W(z') where z' is the point where the antiquark crosses the shock wave. If we substitute 
this quark-antiquark propagator in the original expression for the amplitude (|j) we obtain 

Jd 4 xd 4 ze ipA - x+iq - z (T{j(x + z)j(z)}) A ~ 

/ §K(P±, q±)Tr{U(p ± )W(q ± - p ± )} (11) 
where U(p±) is the Fourier transform of U(x±) and I(p±,q±) is the so-called "impact fac- 



* Because the shock wave is very thin the quark has no time to deviate in the transverse direction. 
Therefore the quark's trajectory inside the shock wave can be approximated by a light-like straight 
line which means that the interaction of the quark with the shock wave will be described by a gauge 
factor ordered along this segment of a straight line. Since there is no field outside the shock- wave 
"wall" one can formally extend the limits of integration in a gauge factor to ±00 which gives the 
operator U 
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tor" which is a function of p±,p±-q±, and photon virtuality fL2 |, |J. Thus, we have repro- 
duced the leading term in the expansion (|T]). (To recognize it, note that lP(x±)U(y±) = 
Pexp {— ig fy dZiUi(z±)\ where the precise form of the path between points x± and y± does 
not matter since this is actually a formula for the gauge link in a pure gauge field Ui(z±)). 
So, in the leading order in perturbation theory we have calculated the integral over fast 
fields explicitly and reduced the remaining integral over slow fields to the matrix element 
of the two- Wilson-line operator, see Fig. f|. It is worth noting that in the next order in 
perturbation theory we will get the contribution to the r.h.s of Eq. (|Tl|) proportional to 
four- Wilson-line operators, in the next to six-line operators and so on. 

Note that formally we have obtained the operators U ordered along the light-like lines. 
Matrix elements of such operators contain divergent longitudinal integrations which reflect 
the fact that light-like gauge factor corresponds to a quark moving with speed of light (i.e., 
with infinite energy). This divergency can be seen already at the one-loop level if one calcu- 
lates the contribution to the matrix element of the two- Wilson-line operator U(x±)W(y±) 
between the "virtual photon states" shown in Fig. B. The reason for this divergency is very 





FIG. 4. A typical Feynman diagram for the 7*7* scattering amplitude (a) and the corresponding 
two- Wilson-line operator (b) 

simple. We have replaced the fast-quark propagator in the "external field" (represented here 
by two gluons coming from the bottom part of the diagram) by the light-like Wilson line. 
In doing so we have assumed that these two gluons are slow, tj p <C tja- However, when we 
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calculate the matrix element of the U(x±)W (y±) formally the integration over the rapidities 
of the gluon r] p is unbounded so our divergency comes from the fast part of the external field 
which really does not belong there. Indeed, if the rapidity of the gluon 7] p is of the order 
of the rapidity of the quark this gluon is a fast one so it will contribute to the coefficient 
function (in front of the operator constructed from the slow fields) rather than to the matrix 
element of the operator. 

This is very similar to the case of usual light-cone expansion for the deep inelastic scat- 
tering (DIS) at moderate x. In that case , we at first expand near the light cone (in inverse 
powers of Q 2 ). The result is that the amplitude of DIS is reduced to matrix elements of 
the light-cone operators which are known as parton densities in the nucleon. These matrix 
elements contain logarithmical divergence in transverse momenta for the same reason as 
above - when expanding around the light cone we assumed that there are no hard quarks 
and gluons inside the proton, but matrix elements of light-cone operators contain formally 
unbounded integrations over k\. It is well known how to proceed in this case: we define the 
renormalized light-cone operators with the transverse momenta k\ > /i 2 cut off and expand 
our T-product of electromagnetic currents in a set these renormalized light-cone operators 
rather than in a set of the original unrenormalized ones (see e.g. ||). After that, the matrix 

2 

elements of these operators (parton densities) contain factors In ^ and the corresponding 
coefficient functions contain ln%. When we calculate the amplitude we add these factors 
together so the dependence on the factorization scale /i cancels and we get the usual DIS 
logarithmical factors In 

Similarly, we need some regularization of the Wilson-line operator which cuts off the fast 
gluons. As demonstrated in 0, It can be done by changing the slope of the supporting line as 
demonstrated in ||. If we wish the longitudinal integration stop at rj — r] , we should order 
our gauge factors U along a line parallel ton = api + ap2, then the coefficient functions in 
front of Wilson-line operators (impact factors) will contain logarithms ~ g 2 In 1/a. Similarly 
to DIS, when we calculate the amplitude, we add the terms ~ g 2 \nl/cr coming from the 



10 



coefficient functions (see Fig. |5|b) to the terms ~ g 2 In -^rji coming from matrix elements 
(see Fig. [|a) so that the dependence on the "rapidity divide" a cancels and we get the usual 
high-energy factors g 2 In -% which are responsible for BFKL pomeron. 






+ 





(a) (b) 

FIG. 5. Decomposition into product of coefficient function and matrix element of the 
two- Wilson-line operator for a typical Feynman diagram. (Double Wilson line corresponds to 
fast-moving gluon) 



III. FACTORIZATION FORMULA FOR HIGH-ENERGY SCATTERING 

In order to understand how this expansion can be generated by the factorization formula 
of Eq. (^|) type we have to rederive the operator expansion in axial gauge A, = with an 
additional condition AJ x = _ oo = (the existence of such a gauge was illustrated in [13| 
by an explicit construction). It is important to note that with with power accuracy (up 
to corrections ~ a) our gauge condition may be replaced by n^A^ = 0. In this gauge the 
coefficient functions are given by Feynman diagrams in the external field 

B t {x) = Uiix^Qix*), B. = B* = (12) 

which is a gauge rotation of our shock wave (it is easy to see that the only nonzero component 
of the field strength tensor F,i(x) = Ui(x±)5(x*) corresponds to shock wave). The Green 
functions in external field (|l2|) can be obtained from a generating functional with a source 
responsible for this external field. Normally, the source for given external field is just 
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J v = D^F^v so in our case the only non- vanishing contribution is J*(B) = D l F*. However, 
we have a problem because the field which we try to create by this source does not decrease 
at infinity. To illustrate the problem, suppose that we use another light-like gauge A* = 
for a calculation of the propagators in the external field (p!2|) . In this case, the only would-be 
nonzero contribution to the source term in the functional integral D l F it A 1f vanishes, and it 
looks like we do not need a source at all to generate the field (This is of course wrong 
since is not a classical solution). What it really means is that the source in this case 
lies entirely at the infinity. Indeed, when we are trying to make an external field A in the 
functional integral by the source J M we need to make a shift A^ — ► A^ + A^ in the functional 
integral 

JVAexp{iS(A) -ifd 4 xJ*(x)A ail (x)} (13) 

after which the linear term D^F^ V A V cancels with our source term J^A^ and the terms 
quadratic in A make the Green functions in the external field A. (Note that the classical 
action S(A) for our external field A = B (IT3) vanishes). However, in order to reduce the 
linear term / d i xF^ v D^A U in the functional integral to the form / d A xD fJ, F fJiU A u (x) we need 
to make an integration by parts, and if the external field does not decrease there will be 
additional surface terms at infinity. In our case we are trying to make the external field 
A = B so the linear term which need to be canceled by the source is 



- fdx.dx*d 2 x±F i .D :¥ A i = jdx^x^F^A^' °° (14) 

x,=— oo 



It comes entirely from the boundaries of integration. If we recall that in our case F,i(x) = 
Ui(x±)5(x*) we can finally rewrite the linear term as 

J d 2 x ± U i (x 1 _){A i {-oop 2 + x±) - A*{oop2 + x±)} (15) 

The source term which we must add to the exponent in the functional integral to cancel 
the linear term after the shift is given by Eq. with the minus sign. Thus, Feynman 
diagrams in the external field (|TJ) in the light-like gauge A* = are generated by the 
functional integral 

12 



I VAexp {iS(A) + % J d 2 x 1 U ai {x i _)[A a i {oop 2 + x±) - A ai (-oop 2 + x ± ))} (16) 

In an arbitrary gauge the source term in the exponent in Eq. ([H]) can be rewritten in the 
form 

p roc 

2i / d 2 x±Tr{U % (x±) / dv[-oop 2 , vp 2 ] Xx F H (vp 2 + x±)[vp 2 , -oop 2 ] Xx } (17) 

J J — oo 

(Hereafter we use the space-saving notation [up 2 ,vp 2 ] x± = [up 2 + x±,vp 2 + x±] and similar 
notation for gauge link ordered along pi). Thus, we have found the generating functional 
for our Feynman diagrams in the external field (13). 

It is instructive to see how the source ([17]) creates the field (|T2|) in perturbation theory. 
To this end, we must calculate the field 

A^x) = J VAA^x) exp [iS(A) + 2iJ d 2 x±Tr{U i (x±) 

[ dv[-<xp 2 ,vp 2 ) x± F*i{vp 2 + x x )[vp 2 ,-<xp 2 } XA }\ (18) 



by expansion of both S (A) and gauge links in the source term ([TtD in powers of g (see Fig. 
0). In the first order one gets 



(a) (b) 




(c) 

FIG. 6. Perturbative diagrams for the classical field JT^ ) 



4 0) (^) = [°° dv [ dz 1 _U ia {z i _){A,{x)F: i {vp 2 + z ± )) (19) 

J —oo J 



where (O) = JVAe tSo O. Now we must choose a proper gauge for our calculation. We 
are trying to create a field (|T2| ) perturbatively and therefore the gauge for our perturbative 
calculation must be compatible with the form ( |T2"D — otherwise, we will end up with the 
gauge rotation of the field B (x). (For example, in Feynman gauge we will get the field A^ 
of the form of the shock wave A{ = A* = 0, A, ~ 5(x*)). It is convenient to choose the 
temporal gauge Ao = [] with the boundary condition A\ t= _ oa = where 

A^x) = f dt'F^{t\x) (20) 

J — oo 

In this gauge we obtain 

A {0 \x) = f dp (g^ - 2 P ^ Pl+P ^ u + ^^ V -l + Ap » Pu \ 



(2tt) 3 V ^ s(a + P + ie) s 2 (a + (5 + ief J 

— — / dz ± e iax ' +ll3x *- l ^- g ^p 2u 5(a-)d i U ia (z ± ) (21) 

aps — p z ± + ieJ 2 

where 5(a|) comes from the / dve lva % . (Note that the form of the singularity ( po ^ e ) which 
follops from Eq. ( p0|) differs from conventional Mandelstam-Leibbrandt prescription V.p.-^). 
Recalling that in terms of Sudakov variables dp = ^dadj3dp± one easily gets that A^ = 
A^ = and 

Aj\x) = 6(x.) J J^y^T J dz^-^-^d^U^) (22) 
which can be written down formally as 

- 0(x»)i$0,-tf J '(*L) = Ui(x ± )9(x*) - 6(x*)^(dl 9ij + d i d j )U j (x ± ) (23) 
o ± o ± 

(in our notations dj_ = —did 1 ). Now, since Ui(x) is a pure gauge field (with respect to 
transverse coordinates) we have dJJj — djUi = i[Ui, Uj] so 



^The gauge .4* = which we used above is too singular for the perturbative calculation. In this 



gauge one must first regulate the external field (12) by, say, replacement Ui6(x*) —>■ Ui0(x*)e ex * 
and let e — > only in the final results. 
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$°\x) = Ui(x ± )9(x*) - e(x*)ig%[U h Uj])(x ± ) (24) 

Thus, we have reproduced the field ( |i~2"D up to the correction of of g. We will demonstrate 
now that this 0(g) correction is canceled by the next-to- leading term in the expansion of 
the exponent of the source term in eq. ([18]). In the next-to- leading order one gets (see Fig. 
|b): 

4° 0*0 =9 JdyJ dz ± dz' ± W a (z ± )U k \z' ± ) (25) 
(A,(x)2Ti{d a Ap(y)[A a (y),A^y)]} JdvF^vpz + z ± ) fdv>F* k (vp 2 + z' ± )) 

It is easy to see that A* = A^ = and 

Af\x)= (26) 
9 JdyJ j0f t e- tp{x - y) ^ (d«[Af\y),AT{y)] + [A^ k (y), d<Af\y) - (i <-> k))) 

Since A k °^ is given by Eq. fl2"4|) this reduces to 

= / (^p^^ &{[UM U k (y)}) + 0(g 2 ) (27) 

which cancels the second term in Eq. fl2"ip. Thus, we obtain 

A(x) = ^(x ± )#(x*) + 0(g 2 ) (28) 

Similarly, one can check that the contributions ~ g 2 coming the diagrams in Fig. |6]c cancel 
the g 2 term in the Eq. fl28|) and so on leading finally to the expression Ui(x±)8(x*) without 
any corrections. 

We have found the generating functional for the diagrams in the external field ([12]) which 
give the coefficient functions in front of our Wilson-line operators C/j. Note that formally 
we obtained the source term with the gauge link ordered along the light-like line which is 
a potentially dangerous situation. Indeed, it it is easy to see that already the first loop 
diagram shown in Fig. [^ is divergent. The reason is that the longitudinal integrals over a p 
are unrestricted from below (if the Wilson line is light-like). However, this is not what we 
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p 



FIG. 7. A typical loop diagram in the external field created by the Wilson-line source ( |l7| ) 

want for the coefficient functions because they should include only the integration over the 
region a p > a (the region a p < a belongs to matrix elements, see the discussion in Sect. 
3). Therefore, we must impose somehow this condition a p > a in our Feynman diagrams 
created by the source (|T7| ) . Fortunately, we already faced similar problem — how to impose 
a condition a p < a on the matrix elements of operators U (see Fig. f|) and we have solved 
that problem by changing the slope of the supporting line. We demonstrated that in order 
to cut the integration over large a > o from matrix elements of Wilson-line operators C/j 
we need to change the slope of these Wilson-line operators to n — op\ + op2- Similarly, if 
we want to cut the integration over small a p < a from the coefficient functions we need to 
order the gauge factors in Eq. ( |T7) ) along (the same) vector n = op\ + ap2 Q. 

Therefore, the final form of the generating functional for the Feynman diagrams (with 

$ Note that the diagram in Fig. ^ is the diagram in Fig. ||b turned upside down. In the Fig. ||b 
diagram we have a restriction a < a. It is easy to see that this also means a restriction (3 > a if one 
chooses to write down the rapidity integrals in terms of /3's rather than a's. Turning the diagram 
upside down amounts to interchange of pa and pb which leads to (i) replacement of the slope of 
Wilson line by ap\ + ap2 and (ii) replacement a <-> (3 in the integrals. Thus, the restriction (3 > a 
imposed by the line collinear to <jp\ + ap2 in diagram in Fig. ||b means the restriction a > a by 
the line || ap\ + ap2 in the Fig. [?] diagram. After renaming a by a we obtain the desired result. 
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a > a cutoff) in the external field (|T3| ) is 

J VAV^exp {iS(A, ^) + iJ d l x ± U ai (x ± )V?{x ± )\ (29) 

where 

Vi(x ± ) = (30) 
S-oo dv[-oon, vn] x± n^F^(vn + x±)[vn, -oon] x± 

and V® = Ti(X a Vi) as usual. For completeness, we have added integration over quark fields 
so S(A, *) is the full QCD action. 

Now we can assemble the different parts of the factorization formula (|j). We have written 
down the generating functional integral for the diagrams with a > o in the external fields 
with a < o and what remains now is to write down the integral over these "external" fields. 
Since this integral is completely independent of ( p9|) we will use a different notation B and 
X for the a < a fields. We have: 

JvAV^e^ A ^j(p A )j(p' A )j(p B )j(p' B ) = (31) 

fDAVi)V^e iS ^j(p A )j(p' A ) PBVxDx 
j(PB)j{p' B )e iSiB ' x} exp {iJd 2 x ± U ai (x ± )V i a (x ± )} 

The operator C/j in an arbitrary gauge is given by the same formula (|3^) as operator Vi with 
the only difference that the gauge links and F 9i are constructed from the fields B^. This is 
our factorization formula (^) in the functional integral representation. 

The functional integrals over A fields give logarithms of the type g 2 In 1 jo while the 
integrals over slow B fields give powers of g 2 \n(os/m 2 ). With logarithmic accuracy, they 
add up to g 2 has/m 2 . However, there will be additional terms ~ g 2 due to mismatch coming 
from the region of integration near the dividing point a ~ a where the details of the cutoff 
in the matrix elements of the operators U and V become important. Therefore, one should 
expect the corrections of order of g 2 to the effective action / dx±U l Vi. Still, the fact that the 
fast quark moves along the straight line has nothing to do with perturbation theory (cf. ref. 
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fl4||); therefore it is natural to expect the non-perturbative generalization of the factorization 
formula (|3T| ) constructed from the same Wilson-line operators C/j and Vi (probably with some 
kind of non-local interactions between them). 

IV. EFFECTIVE ACTION FOR HIGH-ENERGY SCATTERING 

The factorization formula gives us a starting point for a new approach to the analysis 
of the high-energy effective action. Consider another rapidity r]' in the region between 770 
and rjs = lnm 2 /s. If we use the factorization formula (|3l|) once more, this time dividing 
between the rapidities greater and smaller than rj' , we get the expression for the amplitude 
(§) in the form 

iA(s,t) = JvAe tS ^j(p A )j(p' A )j(p B )j(p' B ) (32) 
= JvAe lS ^j(p A )j(p' A ) JvBe iS ^j(p B )j(p' B ) 

J X>Ce* 5(C V J d2x -L vat ( x i-M( x ±)+ i Jd?x ± W ai (x_ L )Y?(x_ L ) 

In this formula the operators Vi (made from A fields) are given by Eq. fl3"0|), the operators 
Ui are also given by Eq. ([30]) but constructed from the C fields instead, and the operators 
Wi (made from C fields) and Yj (made from B fields) are aligned along the direction n' = 
a'pi + a'p2 corresponding to the rapidity 7/ (as usual, In a' /a' = r( where a' = m 2 / sa'): 

Ui(C) x± = dv[-oon, vul^n^F^vn + x x ) [vn, -oon] x± 
Wi(C) x± = J^ 00 dv[-oon',vn'] x± n' IJ 'F IJ ,i(vn' + x±)[vn', -oon'] x± 
Yi(B) x± = J^ 00 dv[-oon',vn%_ L n' tJ 'F tJ ,i(vn' + x±)[vn', -oon% ± 

Thus, we have factorized the functional integral over "old" B fields into the product of two 
integrals over C and "new" B fields. 



§For brevity, we do not display the quark fields. 
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Now, let us integrate over the C fields and write down the result in terms of an effective 
action. Formally, one obtains: 



iA(s,t) = JvAe lS ^j(p A )j(p' A ) JvBe^j( PB )j(p' B y s ^ v > 
where S e s for the rapidity interval between 77 and rf is defined as 



e iS e s(Vi,Yi;fr) _ JJ)C e iS(C) e i J <fx ± V ai {x x )Uf(x ± )+i J d?x ± W ai (x ± )Y?(x ± ) 



(33) 



(34) 



This formula gives a rigorous definition for the effective action for a given interval in rapidity 
(cf. ref. 0). 

Next step would be to perform explicitly the integrations over the longitudinal momenta 
in the r.h.s. of Eq. ([34]) and obtain the answer for the integration over our rapidity region 
(from 77 to rf) in terms of two-dimensional theory in the transverse coordinate space which 
hopefully would give us the unitarization of the BFKL pomeron. At present, it is not known 
how to do this. One can obtain, however, a first few terms in the expansion of effective 
action in powers of Vi and Y{. The easiest way to do this is to expand gauge factors Ui 
and Wi in r.h.s. of Eq. ( |3"4"D in powers of C fields and calculate the relevant perturbative 
diagrams (see Fig. |8|). The first few terms in the effective action at the one- log level Qhave 



w 





(a) (b) (c) (d) 

FIG. 8. Lowest order terms in the perturbative expansion of the effective action. 



the form 



17H: 



** This "one-log" level corresponds to one-loop level for usual Feynman diagrams. Superficially, 
the diagram in Fig. ||d looks like tree diagram in comparison to diagram in Fig. ||c which has 
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(a) (b) (c) (d) 

FIG. 9. Counting of loops for Feynman diagrams (a),(c) and the corresponding Wilson- line 
operators (b),(d) 



S cS = J d 2 xV a \x)Y i a {x)- (35) 
g£a In %(N C J d 2 xd 2 yV%(x) In 2 (x - y) 2 Y^(y) + Jd 2 xd 2 yd 2 x'd 2 y'd 2 z 
V^)V-(y)YUx'm(y') ^ §E$r In (A) 2 In gf^ In gf$) + ... 

where we we use the notation V { j(x) = -^-V^{x) etc. The first term (see Fig. ||a) looks 
like the corresponding term in the factorization formula ([31]) - only the directions of the 
supporting lines are now strongly different fj . The second term shown in Fig. ||c is the 



one loop. However, both of the diagrams in Fig. |8|c and d contain integration over longitudinal 
momenta (and thus the factor m^) so in the longituduinal space the diagram in Fig. [8|d is a 
loop diagram too. It happens because for diagrams with Wilsondine operators the counting of 
number of loops literally corresponds to the counting of the number of loop integrals only for the 
transverse momenta. For the longitudinal variables, the diagrams which look like trees may contain 
logarithmical loop integrations. This property is illustrated in Fig. |9|: the Wilsondine diagram 
shown in Fig. ||b has two loops and the diagram shown in Fig. |9|d is a tree but both of them 
originated from Feynman diagrams shown in Fig. [D|a and c with equal number of loops. To avoid 
confusion, we will use the termin "onedog level" instead of "onedoop level". 
^ Strictly speaking, the contribution coming from the diagram shown in Fig. |8]a has 
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first-order expression for the reggeization of the gluon |l[ and the third term (see Fig. |8|d) 
is the two-reggeon Lipatov's Hamiltonian [p| responsible for BFKL logarithms. 



Let us discuss subsequent terms in the perturbative expansion fl35|). There can be two 
types of the logarithmical contributions. First is the "true" loop contribution coming from 
the diagrams of the Fig.|TD|a type. This diagram is an iteration of the Lipatov's Hamiltonian. 
However, in the same (In ^7) 2 order there is another contribution coming from the diagram 



shown in Fig. 10 b. To treat them separately, we can consider the case when g <^ 1 but the 



(a) (b) 
FIG. 10. Typical perturbative diagrams in the next (ln^-) 2 order. 



sources are strong (~ ^) so gYi ~ gUi ~ 1. In this case, the diagram in Fig.|10|a has the 



order g A Y 2 V 2 (hi^) 2 ~ ( ln ^f while the " tree " Fig-H diagram is ~ g^YfVf i}^^f ~ 
. So, in this approximation the tree diagrams are the most important and should 
be summed up in the first place. As usual, the best way to sum the tree diagrams is given 
by the semiclassical method which will be discussed in next Section. 

However, if we would like to get the result on the one-log level it can be obtained using 
the evolution equations for the Wilson-line operators |J. Note that at this level we have 
only the diagrams of the Fig.O type. These diagrams describe the situation when one of 



the form J d 2 xV at (x)^^-Y a ^ ix) which differs from the first term in r.h.s. of eq. ( |35| ) by 
/ d 2 xV ai (x) ^ (d 2 gij - didj)Y a i (x) . However, it may be demonstrated that this discrepancy (which 
is actually ~ 0(g) for a a pure gauge field Y{) is canceled by the contribution from the diagram 
with three-gluon vertex shown in Fig. ||b just as in the case of perturbative calculation of Ai 
discussed in Sect. 3. 
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the sources is weak and another is still strong (see also refs. pO , ||16|| ). For example, if the 
source Vi is weak (and hence gVi is a valid small parameter) but the source Yi is not weak 
(so that gVi ~ 1 is not a small parameter) one must take into account the diagrams shown 
in Fig. |TT| a and b. The multiple rescatterings in Fig. |TT|a,b describe the motion of the gluon 





(a) 



(b) 





(c) (d) 
FIG. 11. Perturbative diagrams for the effective action in the case of one weak source and one 
strong one. 



emitted by the weak source Vi in the strong external field A{ = Yi8(x*) created by the source 
Yi. These diagrams were calculated in ref. ||. For example, the result of the calculation of 
the diagram in Fig. |TT| a presented in a form of the evolution of the Wilson-line operators C/j 
reads # 



u?( x± ) - y, (*L)- 



WHere Y x = Y(x±) = [oon' + xj_, —oon' + x±] so that Yi(x±) = Y^i^Y x . (Note that we have the 
gauge factors in the gluon (adjoint) representation here). 
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£ In $ J dy ± ^ r (r b %Yjd l Y y ) bc + N C U? 



+ ... 



(36) 



where dots stand for the terms with extra g 2 In ^ factors. This evolution equation means 
that if we integrate over the rapidities rj > rj > rj' Q in the matrix elements of the operator 
Ui we will get the expression (RS) constructed from the operators Y{ with rapidities up 
to rj' Q times factors proportional to g 2 (i]o — r]' Q ) = g 2 \w^. Therefore, the corresponding 
contribution to the effective action at the one-log level takes the form 

(37) 
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8tt 3 



dx ± V i a (x ± )U ai (x ± ) -> / dx ± V?(x x )Y m (x±)- 
]n^fdxj.dy xl ^(i{V i {x ± )Yjd i Y v ) aa - N c V a \x ± )U?(x^ 



where the first term is the lowest-order effective action (= the first term in eq. fl35|)) and 
the second term contains new information. To check this second term, we may expand it in 
powers of the source Yi and it is easy to see that the first nontrivial term in this expansion 
coincides with the gluon-reggeization term in eq. (|35|). 

Apart from the fl3"TD term, there is another contribution to the one-loop evolution equa- 
tions coming from the diagrams in Fig. (Jlljb) M: 



where 



U?{x L )U h Ay A 



4tt 3 



in ± (v? [ I dz ± t^tSt ( Y * Y y + 1 - Y * Y > - Y ' y y 



VtO( XA _) = ^-0(x ± ) - iUi(x x )0(x ± ) 
ox 1 



iO( yA _)Ui(y A 



(38) 



V 



V \ab 



(39) 



are the "covariant derivatives" (in the adjoint representation). The corresponding term in 
effective action has the form 



ig" 



In 



a 



dx ± dy ± (V*Vf) (x x ) / dz 



(x - z) ± -(y-z) 



8vr 3 ~* a' J v * « ' w j — -l ( x _ z y^ y _ ^2 

+ 1 - YlY, - Y z %) ab MV/) (y ± ) 



(40) 



23 



The final form of the one-log effective action for this case is the sum of the expressions 
and 

sSiViM) = J d 2 xV a \x)Y?{x)- (41) 
£ In i Jdx ± dy ±w ^r (^(x^Fj^r - N e V*{x ± )Uf (x ± )) 

where is a weak source and Yj is a strong one. It is clear that if the source V,i is strong 
and Yi is weak as shown in Fig. 0c, d diagrams the effective action S^ J \Vi,Yj) will have 
the similar form with the replacement V «-> Y. 



As we mentioned above, the diagrams in Fig.|10|and Fig. [IT] complete the list of diagrams 
which contribute to the effective action at the one-log level (even if both sources are strong). 
It means that the one-log answer in general case can be guessed by comparison of the 
answers for S^(Vi,Yj) and S^(Vi,Yj) (the simple sum is not enough since some of the 
contributions will be double-counted). Instead of doing that, we will obtain the one-log 
result for two strong sources using the semiclassical method and check that it agrees with 

V. EFFECTIVE ACTION AND COLLISION OF TWO SHOCK WAVES 



The functional integral (|34j) which defines the effective action is the usual QCD func- 
tional integral with two sources corresponding to the two colliding shock waves. Instead of 
calculation of perturbative diagrams (as it was done in previous section) one can use the 
semiclassical approach. This approach is relevant when the coupling constant is relatively 
small but the characteristic fields are large (in other words, when g 2 1 but gVi ~ gYi ~ 1). 
In this case one can calculate the functional integral fl54"|) by expansion around the new sta- 
tionary point corresponding to the classical wave created by the collision of the shock waves. 
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With leading log accuracy, we can replace the vector n by p\ and the vector n' by p 2 - 
Then the functional integral (|53D takes the form: 

e iS eS (Vi,Yi;fr) _ JT>Ae iS< ^ CDiA) ^ S d 2 x ± V a *(x ± )U?(x ± )+i J d 2 x ± W a % a (x ± ) ^g) 

where now 

/oo A /*oo 

dvi^Cvpi+xx), W, a = / d«i^(t^2+ar ± ) (43) 
-oo J —oo 

Hereafter we use the notations 

0{x) = [— oopi + x, x]0(x)[x, —oopi + x] 

0(x) = [—oop2 + x, x]0(x)[x, —oop 2 + x] (44) 

Note that we changed the name for the gluon fields in the integrand from C back to A. 
As usual, the classical equation for the saddle point A in the functional integral ( |42|) is 



5 

6A 



Sqcd + J d 2 x ± V ai (x ± )U?(x ± ) + J d 2 x ± W a % a (x ± 



(45) 



A=A 



To write down them explicitly we need the first variational derivatives of the source terms 
with respect to gauge field. We have: 

5Ui = 5Ai(ocpi + x±) — SAi(— oopi + x±) — /f^ duViSAi(upi + x±) 

6Wi = SAi(oop 2 + x±) - 5Ai(-oop 2 + x ± ) - duV^A^upi + x ± ) (46) 

where 

ViO(x) = diO(x) - i[Ui(x ± ) + A(-oo Pl + x ± ), O(x)} 

V,0(x) = d t O(x) - i\Wi(x ± ) + M-oop 2 + x ± ), 0(x)\ (47) 

Therefore the explicit form of the classical equations ( fi5] ) for the wave created by the collision 
is: 

= (48) 
D^F^ = 5(- s x.)[- s x*pi, -oopi]- s _ L V i V r *(a;i.)[-oopi, fa^pi]^ 
D^F.^ = 5(fx*)[fa;.p2, -oop 2 ] x± 1 ViY t (x ± )[-oop2, ^x.p 2 ] x± 
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Also, as explained in Sect. 3, since our fields do not decrease at infinity there may be 
extra surface linear terms (cf. Eq. (0)) coming from the contributions proportional to 
<L4(±oo) in the r.h.s. of eq. (|46|). The requirement of absence of such terms gives four 
additional equations 

F.i\ Xm=00 = 5(2xJs)Y i (x ± ), F*4*=-oo = S(2x./s)V i (x ± ), (49) 
F.i\ Xt=OQ = 5(2x*/s)[oop 2 , -oop 2 \ Xx Y i (x±)[-oop2,oop 2 \ Xx 
F*i\-. =00 = S(2x,/s)[oop 1 ,-oop 1 ] Xx V i (x±)[-oop 1 ,oop 1 \ Xx 



The two sets of equations ( |4"8| ) and fl4"9] ) define the classical field created by the collision of 
two shock waves 0. 

Unfortunately, it is not clear how to solve these equations. One can start with the trial 
field which is a simple superposition of the two shock waves (|l~2|) 

A (o) = A (o) = q 5 A f) = eMVi + q My . (50) 

and improve it by taking into account the interaction between the shock waves order by 
order |15 |. The parameter of this expansion is the commutator g 2 [Yi,Vk\. Moreover, it 



can be demonstrated that each extra commutator brings a factor In ^ and therefore this 
approach is a sort of leading logarithmic approximation. In the lowest nontrivial order one 
gets: 

/ 'x — z) 2 



x\ + ie 



Af ] =~^j dz x ([Y t (z ± ),V k (z ± )}-z^ k)^-^L\n \l 

= ^-J dz ± ^— ln(-xj + (x- z)\ + ie)[Y k (z ± ),V k (z ± )) 
\bn z J x* + ie " 



where X 1 1 — ~ X ~^ X 9 IS cl> longitudinal part of x 2 . These fields are obtained in the background- 
Feynman gauge. The corresponding expressions for field strength have the form: 



These equations are essentially equivalent to the classical equations describing the collision of 
two heavy nuclei in ref. [jlGf . 
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-Xu 



[x — z) , + ie 



<i) _ 9 



ik 



2tt 2 

(i) _ gs 



F- 



dz i 



-i? + (x — z)^ + 



(x - z)' 



-xjj + (x — z)\ + ie \ x* — ie 



[Y k ,V h 



(K,y fe ]-[y fc ,KD 



57T- 



dz i 



(a; — z)* 



x* + te J 
9ik\Yj,V*] [Y h V k \-[Y k ,Vi\ 



-xl + (x — z)^ 



(52) 



y a;, — ze x, + ie J 

In terms of usual Feynman diagrams (when we expand in powers of source just like in 



previous Section) these expressions come from the diagrams shown in Fig. 12. When we sum 




( a > 0>) (c) (d) 

FIG. 12. Perturbative Feynman diagrams for the field strength ([52]) 



up the three contributions from the diagrams in Fig. 12a,b, and c the three-gluon vertex 



in Fig. |12| a is replaced by the effective Lipatov's vertex and we get fl5"2|) up to the terms 



-^didkY k and ^djd k V k standing in place of Yi and Vj. However, as we have discussed in 



Sect. 3, the difference Yi — jpdid k Y k = g^\Y i: Y k ] (which has an additional power of g) will 
be canceled by the next-order perturbative diagrams of the Fig. |T2] d type. 

Let us now find the effective action. In the trivial order the only non-zero field strength 
components are F,f* = 8(~x*)Yi(x±) and F^ = S(-x,)Vi(x±) so we get the familiar expres- 
sion S<® = / d 2 x 1 V ai Y l a . In the next order one has 

S W = J d'x (-^FP ai F£ )a - l - F ^ a F^ aik + ^F« a F« a ) + 2jd 2 x ± 

J du(TrV l ([-oopi, up 1 ] x± F, i (up 1 + x ± )[up 1 + x ± , -oopi] x . ± ) (1) + 

TtY 1 (-oop2,up2]x ± F*i(up2 + x±)[up 2 , oop 2 \ x _ L )^' ) (53) 



27 



We have seen above that the effective action contains In ^ (see Eq. (|35p). With logarithmic 
accuracy, the r.h.s of Eq. (|53|) reduces to 

S« =-l j d 4 xF^ ai (x)F^ )a (x) 

+ J d 2 x ± 2Tr[Y\ Vi] ([x x , -oop 2 + - [x ± , -oo Pl + xj (1) ) (54) 

The first term contains the integral over d 4 x = ^dx,dx^d 2 x±. In order to separate the 
longitudinal divergencies from the infrared divergencies in the transverse space we will work 
in the d = 2 + 2e transverse dimensions. It is convenient to perform at first the integral over 
x* which is determined by a residue in the point x* = 0. The integration over remaining 
light-cone variable x, factorizes then in the form f£°dx,/x, or J^dx./x.. This integral 
reflects our usual longitudinal logarithmic divergencies which arise from the replacement of 
vectors n and ml in fl34|) by the light-like vectors p\ and p 2 - In the momentum space this 
logarithmical divergency has the form / da /a. It is clear that when a is close to a (or a') we 
can no longer approximate n by p\ (or n! by P2). Therefore, in the leading log approximation 
this divergency should be replaced by In ^: 

[°° 1 f°° 1 f a ' 1 a 

/ dx, — = / da > / da— = hi— (55) 

Jo x, Jo a Jo a a' 



The (first-order) gauge links in the second term in r.h.s. of Eq. ( |54D have the logarithmic 
divergence of the same origin: 

[x ± ,-oo Pl + xJ« = ~l° — [ d 2 x ±T ^^[Y k (z ± ),V k (z ± )} 

S7l z J-oo Xjf J (X — Zfl 

[x ± , -oop 2 + x ± }^ = * f ^ [ d 2 x x . T{€ \ 2e [Y k ( Zl _), V k (z ± )} (56) 

8n 2 J-oo x. J (x — z)f 

which also should be replaced by In ^ . Performing the remaining integration over x± in the 
first term in r.h.s. of Eq. ( |54"|) we obtain the the first-order classical action in the form: 

= (57) 
- g In i / A l( lV(iI(xi)L?(!/i) + ^(xx)^^))^ 



where 
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L a = ^abcyaybj^ = ^yaiybk (gg) 

and 6ik is the totally antisymmetric tensor in two transverse dimensions (ei 2 = 1). One may 
also rewrite this expression in a compact form 

s(1) = f^v/M L ^ L?+LS flW (59) 



A more accurate version of this formula looks like: 
2tt 



In °- J d^(L^L<l + LliY^Y + V^V - ^T b L b 2 + 



_L 

L ^ Y] % Y - Y ^ v ) L ^ lk - L ^ lk ^% Y % - Y ~ v r bL i) + 

0([U,V] 3 ) (60) 

where 

Y(x±) = [oopt, -oopi]a; x , V(x±) = [oop 2 , -ocp 2 } x± (61) 

It is easy to see that in the case of one weak and one strong source this expressions coincides 
with fl40| ) (up to the terms of higher order in weak source which we neglect anyway). 

At d = 2 we have an infrared pole in which must be canceled by the corresponding 
divergency in the trajectory of the reggeized gluon. The gluon reggeization is not a classical 
effect in our approach - rather, it is a quantum correction coming from the loop corresponding 
to the determinant of the operator of second derivative of the action 

5 S 



5A„ 5A U 



Sqcd + J d 2 x ± V ai (xj_)U?(x ± ) + J d 2 x±_W ai Y?{x^ 



(62) 

A=A 



The lowest-order diagrams are shown in Fig. [L3| and the explicit form of the second derivative 
of the Wilson-line operator is: 

SUi = i /f^ du /"^ dv[8Ai(upi + x±), V i 5A i (vp 1 + x±)] 

5Wi = % J^du S^ 00 dv[A i (up 2 + x±),Vi8Ai(up 2 + x_l)] (63) 
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w 



(a) 



(b) 



FIG. 13. Lowest-order diagrams for gluon reggeization. 
Now one easily gets the contribution of the Fig. [13] diagrams in the form: 

S r = (64) 

&\ n ±f d 2 x±d 2y ±{V f(x ± )Y™(y ± ) - V«(x ± )Y™(x ± )) 

A more accurate form of this equation reads: 

S T = (65) 
-^ln- fdVdV F2(1 + e) lmxAY*(xA- 

^-(Y l (x ± ){Y(x ± )Y\y ± ) + V(x ± )V\y ± ) - 1}^))°"} + 0([U, Vf ) 

where O aa =TiO in the gluonic representation. In the case of one strong and one weak 
source it coincides with (]37f) (up to the higher powers of weak source). 

The complete first-order (= one-log) expression for the effective action is the sum of S(°\ 
S'W, and 5 r : 



S ea = I d 2 xV ai (x)Y i a (x)- (66) 
|^ln^/ d 2 xd 2 y{^^(L?(x)L?(y) + L«{x)L b 2 {y){Y^Y y + V}V y - if 

i^ Z-^^ ix)(Y ^ - Y - v)abLliy) - 

L a 2 (x)(YX-Y^Vr b L b M)} 

i(V(a; J .){y( a : J .)yt(y ± ) + V(xj.)y+(y ± ) - 1}^ (y ± )) M } + 0([tf, Vf ) 
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At one weak and one large source it coincides with fl41|). (As we discussed in Sect. 4, the 
new nontrivial terms in the case of two strong sources start from [U, V] 3 In 2 ^). 

As usual, in the case of scattering of white objects the logarithmic infrared divergence 
~ = cancels. For example, for the case of one-pomeron exchange the relevant term in the 
expansion of e lScS has the form: 

- ^ In % J d 2 x L d 2 y L f dam {V«Y m ig tk + V«Y k m - V k a Y™)(x ± ) 

T(e)^ j?dbn ^yfbYnl gik _|_ ^/6i\^mfc ^/^V™) (?/_|_) -|- 

^\n^Jd 2 x ± V^x ± )Y ai (x ± )Jd 2 y x d 2 y' ± (VP(y ± ) - V}(y' ± )) 

ay ^yi } t^ (y bJ (y±)-y bJ (y f ±)) (67) 

It is easy to see that the terms ~ - cancel if we project onto colorless state in t-channel (that 
is, replace V al Vj by j$rzi V cl V?) . It is worth noting that in the two-gluon approximation 
the r.h.s. of the eq. (^) gives the BFKL kernel. 

VI. CONCLUSIONS AND OUTLOOK 

The ultimate goal of this approach is to obtain the explicit expression for the effective 
action in all orders in One possible prospect is that due to the conformal invariance 

of QCD at the tree level our future result for the effective action can be formalized in terms 
of conformal two-dimensional theory in external two-dimentional "gauge fields" Vi and Y^. 

Up to now, we have not used the conformal invariance because it is not obvious how to 
implement it in terms of Wilson-line operators. We can, however, expand Wilson lines back 
to gluons. The conformal properties of (reggeized) gluon amplitudes are well studied now. 
In the coordinate space the BFKL kernel is invariant under Mobius group and therefore 
the eigenfunctions of BFKL kernel are simply powers of coordinates. Moreover, at large 
N c the diagrams with fixed number of reggeized gluons (which form a unitary subset of all 
diagrams) may be described in terms of two-dimensional quantum mechanics of the particles 
with Lipatov's Hamiltonian (|3li|). Due to the property of a holomorphic separability this 
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two-dimensional quantum mechanics reduces to the one-dimesional Heisenberg xxx spin-0 
model |)T5| . (Unfortiunately, the exact solution of this model is not known yet). It is not 
clear which part of this symmetry survives for the full effective action but there is every 
reason to believe that it will simplify the structure of the answer even after reassembling of 
Wilson lines. 

In conclusion I would like to note that the semiclassical approach developed above for the 
small-x processes in perturbative QCD may be modified for studying the heavy-ion collisions. 
As advocated in ref. [pO|l , for the heavy- ion collisions the coupling constant may be relatively 
small due to high density. On the other hand, the fields produced by colliding ions are large 
so that the product gA is not small - which means that the Wilson-line gauge factors V 
and Y are of order of 1. It should be mentioned, however, that in this paper we considered 
the special case of the collision of the two shock waves, namely without any particles in the 
final state. It follows from the usual boundary conditions for Feynman amplitude (^|) which 
we calculate: no outgoing waves at t —>■ oo (and no incoming fields at t — > — oo, but we have 
satisfied this condition by choosing the gauge A\ t ^_ oo = 0). However, people are usually 
interested in the process of particle production during the collision (see e.g. |^T[) since it gives 
the experimental probe of quark-gluon plasma. In this case, our approach must be modified 
for the new boundary conditions — we must solve the classical equations (|48|) with only half 
of the boundary conditions (^) at t — > — oo. The boundary condition at t — ► oo depends 
on the problem under investigation: in the case if we are interested in the wavefuction of 
the system at large times we do not have any boundary conditions at t — > oo but we must 
use the causal (retarded and advanced) Green function instead of the usual Feynman ones. 
( For example, in the expression fl52|) for the field strength we will have the retarded Green 
function 9(x )2ttS(x^ — (x — z)\) instead of the Feynman propagator (— xjj + (x — z)\ + ie)~ 1 ) 
On the contrary, if we calculate the total cross section (cut diagrams) we must calculate the 
double functional integral corresponding to the integration over the "+" fields to the right 



and the "-" fields to the left of the cut (see ref 122]). (This is actually a functional-integral 



formalization of Cutkovsky rules). In this case we may use the usual (Feynman and c.c. 
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Feynman) propagators for each type of the fields. The boundary condition requires that two 
types of the field — the left -side "-" fields and the right-side "+" ones — coincide at t — > oo. 
(This boundary condition is responsible for the S(p 2 )8(po) propagators on the cut). Thus, 
to find the total cross section of the shock-wave collision in the semiclassical approximation 
we must solve the double set of classical equations for "+" and "-" fields with the boundary 
condition that these fields coincide at infinity. The study is in progress. 
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